Bosonization theory for tunneling spectra in smooth edges of Quantum Hall systems 
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We calculate the spectral function of a smooth edge of a quantum Hall system in the lowest Landau 
level by means of a bosonization technique. We obtain a general relationship between the one electron 
spectral function and the dynamical structure factor. The resulting I~V characteristics exhibit, at 
low voltage and temperature, power law scaling, generally different from the one predicted by the 
chiral Luttinger liquid theory, and in good agreement with recent experimental results. 



Edges of Quantum Hall systems have been studied 
mainly as realizations of exotic correlated states of elec- 
trons. For example a sharp edge of a fractional quantum 
Hall fluid is described as a Chiral Luttinger Liquid (CLL) 
[jjj . The main characteristic of a CLL is the lack of single 
particle excitations at the Fermi level: the elementar 
excitations are collective waves. Experimental data 
seem to support this picture at filling factor 1/3. 

In many experiments the potential confining the elec- 
trons is slowly varying on the scale of the magnetic length 
I = {Tic/ eB) 1 / 2 . Under these conditions, the edge re- 
constructs itself j3j. For sufficiently smooth confining 
potential the density profile can be determined (on a 
coarse grained scale > I) by the minimization of a lo- 
cal energy functional (J] . This leads to a picture of the 
edge as a broad compressible region characterized by frac- 
tional, non-uniform occupation of degenerate orbitals at 
the Fermi level. The single particle description is mean- 
ingless for such a state, and a full many-body treatment 
of the excitations is needed. As a first step in that di- 
rection, in a recent paper J5| we have presented a calcu- 
lation of the spectral function of a smooth edge using a 
bosonization technique inspired by previous calculations 
for the uniform electron gas |^| || . The theory contained 
a (physically motivated) cutoff on the number of hydro- 
dynamic modes contributing to the spectral function for 
a given edge width d, and predicted a tunneling exponent 
increasing linearly with d. Similar results have been ob- 
tained by Han and Thouless H[n| using a Green's func- 
tion method at finite temperature. 

In this paper we refine our calculation in such a way 
that short wavelength modes are automatically excluded 
from contributing to the spectral function. Thus, the 
present theory is parameter free. We also extend our 
calculations to finite temperature fll| and obtain the low 
voltage/temperature scaling of the spectral function and 
of the tunneling current. Our numerical results, in the 
sharp edge limit, are found to be in good agreement with 
a recent experiment by Chang et al. jfj] . 

Let us begin by writing down the microscopic Hamil- 
tonian within the LLL in terms of density fluctuations 
relative to the equilibrium density profile po(y). 
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where p(f) = 8p{r) + po(y) is the density operator pro- 
jected in the LLL. The hamiltonian ([!]) is diagonal in 
terms of the operators |Bj 



dxe~ ikx / dyf nk {y)Sp(x,y), (2) 
y'kl^pL Jo Jo 

where f n k (y) are the solutions of the equation 



K (k\y-y'\)f nk (y')^-dy' = -?-/„* (j/), (3) 
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and Ko(y) is the modified Bessel function. They form 
a complete set and satisfy the orthonormality condi- 

tion f« f nk (y)f mk (y)£!M. 



-dy = S nm , and vanish out- 
side the interval [0, d\. Equation (|J) is equivalent to 
the eigenvalue problem obtained by Aleiner and Glaz- 
man []13[ from their hydrodynamical treatment of the 
collective modes of the edge. The eigenfunction of the 
n-th branch has n nodes in the y direction and energy 
Unk = kve 2 / X nk Tr, where v = 2irl 2 p is the usual filling 
factor in the bulk. It has been shown in Ref. || that the 
operators b nk satisfy boson commutation relations in the 
long wavelength limit (hi <C 1 and n <C d/l). 

The electronic spectral function can be obtained within 
the independent boson model (IBM) p|,|l4|], which de- 
scribes a single localized electron electrostatically cou- 
pled to density fluctuations. The IBM hamiltonian is 



Hibm= ^2 huj ^bl k b nk + 

nk>0 

+^(r)m ]T M nk (y)[bl k e ik * + b nk e- ik % (4) 
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where ip' (r) is the field operator that creates an electron 
in the LLL coherent state (gaussian) orbital centered at 
r, which is coupled to the bosons by the matrix element 

M nk (y) = 2 -f f dy'e-^l l \~^l^ nk {y')M. 
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(5) 

The form of the matrix elements differs from the one 
used in Ref. || by the presence of the gaussian convo- 
lution. The main consequence is that M nk {y) is found 



f 



to vanish exponentially for large n and/or k, and there- 
fore short wavelength modes are automatically excluded. 
The hamiltonian (|J) can be solved by standard methods 
fbl , within the one-electron Hilbcrt space. The fermionic 
Green's function G > (y; t) = —i(ip(f, t)ip (r, 0)) is propor- 
tional to 



M nk(y) 



v7l&>0 



(N nk + l)(e" 



1) 



+ N nk {e 1 ' 



1) 



(6) 



where N nk is the Bose occupation factor. The local den- 
sity of states is given by the spectral function A > (y,co), 
which is defined as the Fourier transform of iG> (y, t) /2ir. 
From cq. (||) it can be easily shown that A > (y, ui) 
satisfies the integral equation 



wA>(y,u;) 



where 
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A>(y,a;-fi)dn, (7) 
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[5(0 - Wnfc ) + <5(n + Wnfc)] . (8) 



In general, g y (£l) can be obtained from the imaginary 
part of the density-density response function x(r, f , f2), 



ff „(fi) =- /" d 2 r'd 2 r"v(r-r') U (r-f") 



Imx(ry,Q) 
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where u(r) has the Fourier transform w(fc) = 
(2we 2 /k) exp (— (kl) 2 /A). This important connection be- 
tweeen the spectral function and the dynamical structure 
factor appears to be a general property of quasiclassical 
collective states in the LLL. In the case of a uniform 
electron liquid in the LLL equation (|^) was first written 
down by Haussmann JtJ . 

If the function g y (fl) has a finite limit for — > 0, 
at sufficiently small to and T (u>,T <C ojq) the spectral 
function obeys the scaling relation 

A > ( SlW )«T*C«- 1 /( ] ^ I j»(fl)) (10) 

where [[01 f(x,a) has a finite limit for x — > and behaves 
as x Q_ Mbr x ^> 1; in particular fix, 2) = x/(l — e _a: ). 

The tunneling current /(V) for vertical tunneling ex- 
periments between a smooth edge and a metallic contact 
satisfies a similar scaling relation J16|, 



I{V) oc T g ^ Q) h 
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where 



h(x,a) = / f(y,a) 



1 
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is linear at small x and proportional to x a at large x. 
Therefore Ohm's law will be satisfied for V <C fcsT, while 
for V larger than ksT but still small on the electronic 
energy scale uj one gets a power law I(V) ~ T/ 9 "^ . The 
normal Fermi-liquid behaviour is recovered for g y (0) = 1 
(h(x, 1) = a; at any x). 

The calculation of the exponent g y (0) is easily per- 
formed. Neglecting the weak nonlinearity of the n = 
mode we obtain a — g v (0) — J2 n 0n(y), where 
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My) = -f n0 (y) e 



-n 2 /2m 2 



(13) 



and m = d^/2/tir. The gaussian factor accounts for the 
convolution with the gaussian wavefunction in equation 

We observe that independently of the shape of the den- 
sity profile /3o(y) — with negligible corrections aris- 
ing from the weak nonlinearity of the dispersion of the 
ri = mode. In the sharp edge limit, when only one 
branch of edge waves exists, we obtain A > (uj) ~ uj v I v ~ x . 
For v = 1/3 this agrees with both CLL theory and ex- 
periment 0. The exponent g y (Q) — d/lv increases lin- 
early with d and therefore in the limit d — > oo (limit of 
infinitely smooth edge) the tunneling density of states 
vanishes at low energy faster than any power law. This 
is consistent with previous results for the uniform system 
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FIG. 1. Electronic spectral function A > (u) as a function 
of lu/u>o, where ujq — ve 2 /nl, for edges of a v = 1 QH system 
with m — 1, 2, 3, 4, 6 and 12, at an intermediate position 
y = 0.6d. Logarithmic corrections to the edge magnetoplas- 
mon dispersion are here neglected. 
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FIG. 2. Tunneling current I(V) between a metallic contact 
and a sharp edge of a v — 1/2 QH fluid at various tempera- 
tures, on a bilogarithmic scale. Full curves give the present 
result, points are the experimental data by Chang et al. Jl2| . 
The inset shows the corresponding small-oj spectral function 
A > (uj), on a linear scale, compared with the low and high 
frequency asymptotic behaviours (dashed lines). 

In Figure [I] we present our numerical results for the 
full electronic spectral function, calculated from eq. (Q) 
for different edge widths d. Notice that the calcula- 
tion is free of adjustable parameters. The eigenfunc- 
tions f n k(y) and eigenfrequencies uj n k have been calcu- 
lated solving numerically equation (^) at finite k in the 

lOl. This form of 



iasn^ 



special case po(y) = p 
Po(y) allows an analytical solution at small k, which gives 
My) = \Tl (2^) (2 - <We- 2/2m2 ~ e- 2 / 2 ™ 2 M 
where T n (y) is the n-th Chebyschev polynomial. 

Figure [2| compares our low-frequency results with the 
experimental data by Chang et al. |12| for tunneling 
between a metallic contact and a nominally atomically 
sharp edge of a v — 1/2 QH system. The current has 
been obtained from equations (11- 12]) assuming that only 
the n = charged mode contributes to the response, i.e. 
with an exponent a — 2. Both the Ohmic behaviour 
at low voltages and the non-Fermi liquid behaviour at 
larger voltages are in good agreement with experiment. 
In contrast to this, a recent composite fermion theory by 
Shytov et al. |18| predicts a = 3 in this case. 

In conclusion we have obtained via a bosonization ap- 
proach a general relation between the one-electron spec- 
tral function and the dynamical structure factor. This 
allows a parameter-free calculation of the spectral func- 
tion which interpolates between the sharp edge limit and 
the uniform electron gas and is in good agreement with 
experiments on sharp edges both at compressible and in- 



compressible filling factors. 
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